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Abstract: We propose using the extra dimension separating the domain walls carrying 
lattice quarks of opposite handedness to gradually filter out the ultraviolet fluctuations 
of the gauge fields that are felt by the fermionic excitations living in the bulk. This 
generalization of the homogeneous domain wall construction has some theoretical features 
that seem nontrivial. 
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1. Introduction. 

To explain our basic idea it is best to start using continuum notation and language. 

In Callan and Harvey's original domain wall construction one had fermions propagating 
in a full five dimensional gauge background. This situation was maintained in Kaplan's 
original proposal, but was early on discarded in favor of a five dimensional gauge field that 
is restricted to be independent of the fifth dimension and whose component in the fifth 
dimension is zero. 

Here we focus only on the vector-like case, with two domain walls, separated by a 
fifth dimension of the topology of a circle. The domain walls are at diametrically opposed 
points on the circle. On the lattice, one can take the fermion mass to infinity on one of the 
semicircles, and effectively reduce the circle to a segment, with a domain wall at each end. 
The Weyl partners that make up a Dirac fermion live separately on the walls, with small 
leakage into the "bulk" , the interior of the fifth dimension. For a finite separation there is 
an exponentially small direct interaction between the partners, which can be thought of as 
an effective vertex of the structure of an exponentially small mass term. Sometimes, in the 
lattice field theory literature, such mass terms are referred to as a "residual mass" . Only at 
infinite separation does one get exactly massless Dirac particles, with the associated exact 
chiral symmetry. The latter holds since one can independently rotate the partners due to 
the infinite separation between the walls. 

Here we shall re-introduce a dependence of the four components on the fifth dimension. 
For the time being, we maintain the fifth component at zero. However, the five dimensional 
gauge field is uniquely defined in terms of a four dimensional background; the dependence 
on the fifth dimension is chosen by hand, in a way designed to improve the speed of 
convergence to an infinite separation limit. 

For our construction we need to introduce the concept of a "UV" filter. This is a term 
making its appearance in the lattice field theory literature, but its meaning varies slightly, 
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depending on context. To be specific, we first define the term in the continuum. A "UV" 
filter T will be an operation that gets as input some four dimensional Euclidean gauge 
field A and produces a new four dimensional Euclidean gauge field T{A). We require that 
T and the operation of a four dimensional gauge transformations of A commute. J 7 (A) is 
smoother than A in that its the gauge invariant content fluctuates less than that of A. As 
an example, we define a filter using APE smearing: 

F T0 (A) = A(r ), d T A^r) = d A ^S{A)\ A=A{r) , A(0) = A (1.1) 

Here S is some local four dimensional gauge invariant action, like the ordinary Euclidean 
Yang Mills action possibly including also higher derivative terms, ro is a quantity of 
dimensions length (for separations along the four physical dimension) squared, ^/tq is a 
relatively short distance on the scale of four dimensional physics. Linearizing the above 
equation is required to produce, for r > a suppression for higher momentum modes in 
A(t) in Feynman gauge. 

We want to use a UV filter on the gauge fields seen by most of the fermions (bulk) to 
hasten the convergence to the infinite wall separation limit. However, we do not want the 
UV filter to affect too much the physical fermions, which live on the walls, because, for 
example, if we add a mass to make the physical quarks heavy, we still want the charmonium 
spectrum to come out right, and this requires the Coulomb potential to be well represented 
even at short distances, making UV filtering undesiarble. 

This leads us to introduce a profile in the five dimensional gauge field. We label the 
segment connecting the domain walls by the variable s, where < s < S and introduce 
the profile function r(s) which obeys r(s) = t(S — s) and r(0) = 0. r(s) increased from 
gradually to some value ro, stays there for a stretch of s given by I, until it reaches 
the midpoint s = S/2, after which it reverses its behavior in accordance with the above 
mentioned constraint under reflection. This constraint is needed to produce an extra 
symmetry, a parity operation connecting the fermions bound to the walls. When I = oo 
the walls are no longer coupled, and we have exact chiral symmetry. 

As is well known, it is useful to view s as an Euclidean auxiliary time, with an associated 
s-dependent Hamiltonian H(t(s)). Obviously, the s-independent case is easier to analyze. 
What we gain here is that the gap, g, separating negative and positive energies in H(s) 
increases in the "flat" region of length I. The approach to infinite I is governed essentially by 
e - !^, so it becomes obvious why we want g to be large. The impact of the UV fluctuations 
contained in A(t) is such that \g\ increases with r. 

There can be other choices of the filter and profile. Also, the mass term can have a 
profile, like in the original constructions. What was described above is close to what we 
actually implemented on the computer, in a lattice version of the above. 

2. Definitions. 

We shall label the four dimensional slices by s = 1,2, the physical fermions live on 
slices 1 and S and the associated five dimensional wave functions are localized in the fifth 
coordinate s, to the vicinities of these walls. 
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Following 1T|, we write the kernel for the five dimensional fermionic action in the 



following form: 
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The matrix D is of size 2k x 2k where the entries are q x q blocks and k = S and 
q = 2NL A , where the gauge group is taken as SU(N). The factor 2N counts spinorial and 
gauge group indices. 

The matrices B s and C s are dependent on the gauge background defined by the collec- 
tion of link matrices U^(x). These matrices are of dimension N x N. \i labels the positive 
4 directions on a hypercubic lattice and U^(x) is the unitary matrix associated with a link 
that points from the site x in the /2-direction on slice s 

(C s ) X ai,yf3j = \ X^=l v'/l' 3 [3y,x+fl(Uji(x))ij — 5 X) y+fr(Uji {y))ij} = Y^=l a ^ XWjj)xi,yj 
(B s o)xai,y(3j = \^afi Y^i=l[^xy^ij ~ 8y,x+fi(U*(x))ij — S x ,y+fi,(Up (y))ij] 
{Bs)xai,yf3j = (.Bso)xai,y(3j ^Q^xai,y/3j 

(2.2) 



The indices a,j3 label spinor indices in the range 1 to 2. The indices i,j label color in 
the range 1 to 2. The Euclidean 4x4 Dirac matrices 7^ are taken in the Weyl basis where 
their form is 

7, = ("» (2.3) 







As long as Mq > the matrices B s are positive definite due to the unitarity of the 
link variables. To make almost massless quarks on the lattice one also wants |Mq| < 1. 

It is convenient to introduce the 2q x 2q matrix T5 representing the regular 75 matrix 
on spinorial indices and unit action on all other indices. In terms of q x q blocks we have 
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3. The effective four dimensional fermion action. 

Following the method of p] we get 



det-D = (-) qk ([[detB s )det 



8=1 
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Define 



T, = e 



-H s 



^ S W S C s ^-C s + B s 



(3.2) 



For any s, detT s = 1. 

T; is a symmetric product of T~ l factors. By definition, gauge fields and mass param- 
eters Mq labeled by s and S — s are identical. 



T, = R^T~ il+2) R 



(3.3) 



The R, i?' factors are complex matrices (neither unitary, nor hermitian) representing 
"ramps" and are given by: 



R 



rp—lrp—1 rp — 1 
l r l r _ 1 ..l 1 , 



ATI — 1 rri — 1 m — 1 

1 S 1 S-V 1 S-r 



(3.4) 



In between the "ramps" we have a uniform "plateau" with identical transfer matrix 
factors. The total number of slices is k = I + 2r 



det D = R 9fc (n det B a ) det [l + Tr 1 ] det 



i + r. 



l-T; 



(3.5) 



The same derivation for Pauli Villars (PV) fields should give a PV determinant made 
out of the first two factors in the formula above. Dividing the two expressions, leaves us 
with the last factor as representing the almost massless Dirac fermion. 

Thus, the operator 



Ei 



1 + T, 



(3.6) 



carries the I dependence and its spectrum governs the approach to the chiral limit I = oo. 

Unlike in the homogeneous case, we are so far unable to prove that E\ indeed has a 
limit, -Eqo, as I tends to infinity. If we assume that such a limit exits, the rate of convergence 
would be controlled by the eigenvalue of T, that is closest to unity as I becomes very large. 

Physically, using effective Lagrangian intuition, one would expect E'oo to exist for 
most four dimensional gauge fields. Simulations with the five dimensional action could 
practically resolve this question. 



4. How to test the method. 

Before proceeding to a five dimensional implementation one would like to find by numerical 
methods the lowest or few lowest eigenvalues of the matrix T; + T, _1 but this is tough 
because the condition number of T/ rapidly becomes too large as the number of slices 
increases. A possible trick to get around this goes as follows: 
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It is easy to check that the spectrum of 
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is given by A n) fc = p n e s e s , with A; = 0,1,2,... ,5 — 1 where A n = p n z are the 

eigenvalues of = A\....A S . In our application, II5 is hermitian and positive, so 4> n = 0. 

1/5 

From the relation |A n fe| = p n we see that we would get directly the quantity we expect 
to have a finite limit as I — > 00. Writing p n = e SK ", with a weak 5-dependence in n n 
suppressed, we can find the dominating K n using a routine based on a restarted Arnoldi 
procedure looking for the eigenvalues of Bs — 1 of smallest magnitude Q . 

In some sense, it might be surprising that this relatively simple generalization of the 
homogeneous domain wall setup is substantially more difficult to analyze theoretically, with 
or without numerical tools. Since would provide a new variant of an overlap operator, 
the issue might be of interest separately from the question of increased numerical efficiency 
of QCD domain wall simulations. 

5. Further directions. 

One could adapt to the lattice a definition of A(s) based on simple interpolation. This 
would allow for an adiabatic limit, with an almost homogeneous extra dimension where 
theoretical issues would be under better control, setting aside the question how far one can 
deviate from this limit. 
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